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The Kadanoff-Baym equations (KBE) are usually derived under the assumption 
of the weakening of initial correlations (Bogolyubov's condition) and, therefore, 
fail to correctly describe the short time behavior. We demonstrate that this as- 
sumption is not necessary. Using functional derivatives techniques, we present a 
straightforward generalization of the KBE which allows to include arbitrary initial 
correlations and which is more general than previous derivations. As a result, an 
additional collision integral is obtained which is being damped out after a few col- 
lisions. Our results are complemented with numerical investigations showing the 
effect of initial correlations. 



1 Introduction and basic equations 

Nonequilibrium properties of many-particle systems are traditionally described 
by kinetic equations of the Boltzmann type. Despite their fundamental char- 
acter, these equations have well-known principal shortcomings, e.g. (i) the 
short-time behavior (t < T corr ) cannot be described correctly, (ii) the kinetic 
or the quasiparticle energy is conserved only instead of the total (sum of ki- 
netic and potential) energy, (iii) no bound states are contained, and (iv) the 
equations yield the equilibrium distribution and thermodynamics for ideal par- 
ticles. 

Since the 1950ies, there were many efforts to overcome these deficien- 
cies. In the time-diagonal formalism, Prigogine and Resiboisu, ZwanzigB, and 
Zubarevu derived exact kinetic equations which (i) are nonlocal with respect 
to the time (non-Markovian), (ii) contain an initial correlation contribution, 
and (iii) yield the correct conservation laws for nonideal systems. On the 
other hand, Kadanoff and Baymcl, and KeldyshH developed equations for the 
two-time correlation functions. These well-known Kadanoff-Baym equations 
(KBE) have a non-Markovian character, too. The structure of approximations 
is determined by a single central quantity, the self-energy, and conservation 
laws for nonideal systems are fulfilled for any conserving approximation for 
the self-energyQ. However, the original KBE contain no contribution of initial 
correlations and so do not reach the degree of generality e.g. of the equation 
of Prigogine with respect to this feature. 

In order to solve this problem, several methods have been used, e.g. an- 
alytical continuation to real times which allows to include equilibrium initial 
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correlations 000 and perturbation theory with initia]_correlations0'0@. A 
convincing solution has been presented by Danielewicz u. He developed a per- 
turbation theory for a general initial state and derived generalized KBE which 
take into account arbitrary initial correlations. A straightforward and very 
intuitive method which is not based on perturbation theory Ej'Lj, uses the 
equations of motion for the Green's functions. 

Starting point of this non-perturbative method is the Martin-Schwinger 
hierarchyEil, a system of coupled equations for the s-particle Green's functions 
defined on the Keldysh-Schwinger contour C by 

g x ... s {l...s, l'...s') = (J\ S (T C [*(l)...vl/( s )*+( s ')...* + (l')] ) , (1) 

where 'J is the field operator, Tc the time ordering operator on the contour, 
and (...) denotes averaging over the density operator p. The first hierarchy 
equation is the equation of motion for the single-particle Green's function 



dl{ 3 r\l,l)-[/(l,l)} 5l (l,l';[/) 

- 6(1 - 1') ± i j d2 V{\ - 2). 9l2 (12, l'2+), (2) 
with U being an external potential and g\ the inverse Green's function 

•5"'M-('BT + S) ,(1 - 1 ' ) - (3) 

Eq. (||) is not a closed equation for gi(l, 1'), because the interaction leads to 
a coupling of the one-particle to the two-particle function and so on. 

To find a closed equation for the one-particle Green's function, one has to 
consider three main steps: 

I The hierarchy can be decoupled formally by introduction of the self- 
energy, taking into account, that g\2 can be derived from g\ by means of 
functional derivation, 

J taS(l,I)ffi(l,l , )=±i / d2V{l- 2) 9 i 2 (12, 1'2+) 
c 

= ±ij d2V(l - 2) j± MWp + gi(1; V)gi% 2+) | . (4) 
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II An important consistency criterion for E is the correct limit for t — > to. 
From Eq. (Q), it is obvious that 

lim / tfi 2(1,1)^1(1,1') =±i [ dr 2 V{r l - r 2 ) [c(nr 2 , r' l5 r 2 ; t ) 
ti=t[— >t J J 
c 

+ gi(rir' 1 ,ta)gi(r 2 r' 2 ,to) ±ffi(rir^,to)sfi(r 2 ri,to)] , (5) 

where c(rir2, r^r^; to) is the initial binary correlation. Therefore, £ must 
contain, in addition to the Hartree-Fock contribution, another time-local 
part S m , 

S(l, 1') = E HF (1, 1') + S c (l, 1') + E m (l, 1'), 
E in (l,lO = E in (l,r' 1 *o)«(ti-to), (6) 
and analogous for the adjoint equation, 

S(l, 1') = £ HF (1, 1') + S c (l, 1') + E m (l, 1'), 
£ in (l, 1') = S m (nt , l')S(h - t ). (7) 

III In order to make £ unique, initial (or boundary) conditions are necessary. 
For example, in thermodynamic equilibrium, the solution is uniquely 
fixed by the Kubo-Martin-Schwinger condition. Furthermore, in the 
nonequilibrium case, we can use the familiar Bogolyubov condition of 
weakening of initial correlations, 

lim 5l2 (12,l'2')| t0 = [ffi(l,l')5i(2,2')±.9i(l,2')5i(2,l')]| t0 . (8) 

to — oo 

In this case, E iri = 0, and we recover the orig inal KBE0. 
The most general initial condition has the form 

lim g 12 (l2, 1'2') = [<?i(l, l')ffi(2, 2') ± 9l (l, 2') 9l (2, l')]| to 

t\— 12— 1 ± — t 2 ^to 

+c(r 1 r 2 ,r' 1 r 2 ;to). (9) 
2 Self energy and initial correlations 

In order to determine the self-energy, we start from Eq. ([!]) . With this defini- 
tion we have obtained a formally closed equation for the one-particle Green's 
function on the contour which may be cast into the form 



dl 



(1, 1) - U(l, 1) - E(l, !)} gx (1, 1', U) = 5(1 - l'). (10) 
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This equation is a compact notation of the Kadanoff-Baym equations and 
is sometimes called Dyson-Schwinger equation. The same procedure for the 
adjoint of equation (g) leads to the adjoint Dyson equation with the self-energy 
£. Eqs.(|,@) show that E = £ for all times t, t' > to- 
Further analysis of Eq. (^) requires to evaluateJhe functional derivative 
5g/SU, for which a simple procedure has been givenEl. 

Starting from the Dyson equation, which, for t,t' > to, can be written in 
the form 

J di ffl - i (i,i) ffl (i ) i')= s(i -n (ii) 

c 

where we introduced the inverse Green's function 

^(1, 1') = .9?" (1, 1') - U(l, 1') - E(l, 1'), (12) 
one obtains easily by functional differentiation 

rfI ^rHM) (I V) = _ f di 1(1 ijW. (13) 

SU(2',2) y v ' ' J yi v ' ' 6U(2',2) v ; 

c c 

Using (|l2|), the general solution of this equation and its adjoint is found to be 



MM') 



/• = o E(l,l) + S m (l,l) 

6U(2',2) =^) 9lM+ J d UT gi{ l,l) L i ffl (T,l') 

c 

±C(12,1'2') = ±L(12, 1'2'), (14) 
where C is an arbitrary function which obeys the homogeneous equation, i.e. 

J dl 9 i\l, 1X7(12, 1'2')=0. (15) 
c 

There are three similar conditions, one following from the crossing symmetry 
(1 <-> 2) and two from the adjoint Dyson equation. 

Let us now analyze the physical and mathematical consequences of the 
function C(12, 1'2'). To this end, we consider Eq. (14) in the limit t,t' — » <o- 
In this case, the integral over the contour vanishes, and it directly follows 

L(r 1 ,r 2 ,r' 1 ,r' 2 ,t ) = c(n, r 2 , r\ , r' 2 , t ) ± g 1 (r 1 , r' 2 , t )g 1 (r 2 , r[, t ). (16) 

Hence, the function c(to) is to be identified with initial binary correlations. 
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In order to explore the temporal evolution of the function C(12, 1/2'), we 
solve Eq. (15) (and the three analogous relations) together with the initial 
condition 

a(12 1 l'2')| il=t2 = tl=ti = 4o =C(t ). (17) 

The result is 

C(12,1'2') = J dldMJ 5l (l^) 5l (2,2)c(T2,r1) 5l (f,2 , )<?i(T,l'),(18) 
c 

with c(12, 12) = c(fito,Mo,^i*o,^2*o) 

xS(ti - t )6(t 2 - to)5fa - t )S(h - to). (19) 

Let us now come back to the self-energy. Introducing Eqs. ( |l4| ) and (18) 
into (Q), we find the following functional equation for S: 

^[E(l,l')+S in (l,l')] 



£(!,!') = ±ij d2V(l - 2) i ±J dl.gi(l,l 



5U(2+,2) 

+<5(l-l')5i(2,2+)± ( 5(2-l').gi(l,2 + ) 



dld2d2 9l (l, i) 9l {2, 2) c(12, 1'2) 3l (2, 2 




An analogous equation follows readily for S. With Eq. (|20|), the self-energy 
is given as a functional of the interaction, the initial correlations and the one- 
particle Green's function. From the definition of c , Eq. (p^), it is obvious that 
the last contribution on the r.h.s. is local in time with a <5-type singularity at 
t = t' . Additional terms of this structure arise from the functional derivative. 
A further important property of the self-energy follows from comparing S, Eq. 
( ^p| ) with the corresponding expression for E. One verifies that S = £ for all 
times t, t' > to, which means, in particular, that for these times, a well defined 
inverse Green's function does exist. 

Eq. (|2Cj) is well suited to come to approximations for the self-energy. By 
iteration, a perturbation series for £ in terms of g, V and C can be derived 
which begins with 

E x (l, 1') - ±iS(l - 1') J d2 V(l - 2) 9l (2, 2+) 

±i J d2V(l-2) J dld2d% gi (1, l)ffi(2,2) c(12,l'5) g x % 2+) 
c 

+ exchange. (21) 
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It is instructive, to represent formula (21) in terms of Feynman diagrams: 




where, in contrast to conventional diagram techniques, we have introduced the 
initial correlation as a new basic element, drawn as a shaded rectangle. Second 
order contributions are evaluated straightforwardly too, with the result 




+ ... + exchange. (23) 

The analysis of the iteration scheme allows us to conclude that all contributions 
to the self-energy (all diagrams) fall into two classes: (I) the terms E ffF und 
S c which begin and end with a potential and (II), E m - those which begin 
with a potential, but end with an initial correlation. This conclusion verifies 
the structure of the self-energy, Eqs. (||,0), as we discussed in the previous 
section. The initial correlation part of the self-energy in fact turns out to be 
temporally local (similar to the Hartree-Fock term) which is nonzero only if t[ 
(or, in the adjoint case, ii) is equal to the initial time tg. The same result was 
obtained by Danielewicz based on his perturbation theory for general initial 
states 0. 

If one considers the first two iterations for the self-energy, Eqs. ( p2| , p2| ) , 
more in detail, it becomes evident that, in the initial correlation contribution, 
in front of the function c, appear just the ladder terms which lead to the 
buildup of the two-particle Green's function. Thus, obviously, the iteration 
"upgrades" the product of retarded one-particle propagators in the function 
C to a full two-particle propagator, in the respective order, i.e. £ m is of the 
form 

E'"(ll') = ±i J d2V i l - 2 ) / dridr 2 dfidf 2 

x ff f 2 (12, fit , ratoMfiio, f 2 t a , I'^ioWC^o, 2 + )<5(t' 1 - t ).(24) 

Notice especially, that this renormalization occurs not symmetrically, i.e. in 
the adjoint of Eq. J24|), = [S m ] appears the adjoint propagator g^ 2 . 
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3 Generalized Kadanoff Baym equations 



Let us now come back to the Kadanoff-Baym equations. In order to discuss 
the influence of initial correlations, we insert expression (||) into Eq. (^), 



dl 



(11)-E/(11) 



(11) fll(ll') 



= 5(1 - 1') + J dl [S c (ll) + E m (ll)] (25) 
c 

from which we obtain the KBE for the correlation functions, if we restrict the 
time arguments to opposite branches of the contour: 



^5^(11)^(11') 



dl T, R (11) g< (!!') + dl E*(ll) + E in (ll) gf (11'). (26) 



In contrast to the original KBE, there are two important new properties which 
have to be underlined here: The equations (|2^) are valid for an arbitrary 
initial time point to, and they explicitly contain the influence of arbitrary initial 
correlations in the additional self-energy term E m . The analytical properties 
of the retarded and advanced Green's functions give rise to a damping 712 
leading to a decay of the initial correlation term after a time of the order t ~ 
1 /712 ~ T cor . Thus, there is no need at all to postulate Bogolyubov's weakening 
condition; for t > r cor , the generalized Kadanoff-Baym equations switch from 
the initial regime into the kinetic, or Bogolyubov regime, "automatically". 



4 Non Markovian kinetic equation in binary collision approxima- 
tion 

In order to compare our results with previous work in the field of kinetic theory, 
we will derive a non-Mar kovian kinetic equation from the generalized KBE. 
To get such an equation for the single-time Wigner distribution, the following 
steps are necessary: 

• The KBE ( p6| ) and its adjoint have to be taken in the equal time limit 
and substracted from each other to yield the time-diagonal KBE. 

• Since this equation is not yet a closed equation for the Wigner function, 
the reconstruction problem = g<[f] has to be solved. Wejjse the 
generalized Kadanoff-Baym ansatz proposed by Lipavsky et al. 
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• An approximation for the self-energy has to be introduced. Here, we 
take S in binary collision (T-Matrix) approximation. 

prforming these three steps, we get the non-Markovian Boltzmann equation 

T 



V fl /i(T) = I 1L -(T) + Tr dtdt dt 



t 



x {T R (T, t) U A {t, t) T A (i,i) U A {i,T) [F < (i)F > (t)-F > {t)F < (t)~ 

+ T R (T,1) U R (t,l) T A (t,t) U A (t,T) [F < (t)F > (t)-F > (F)F < (t)]} 
+adjoint (27) 

with 

F< = /!< / 2 < =fih, F>= f> /> = (1 ± h) (1 ± h) 

(-) R/A <-) , „ 

U {tt') = Q[±{t-t')]e lE ^ t - t ^ 

The initial correlation contribution reads 

I IC {t) =nTrl / diT R {t,t) U° R (t,t ) c(t ) U° A (t ,t) 



dtdtdt g° 12 (t,t)T R {t,t)U (t,to)c(t )U° (t , t) T A (t ,t) 
+adjoint (28) 

Here, T R / A (t, i') is the off-shell two-particle scattering matrix given by a gen- 
eralized Lippmann-Schwinger equation. 



5 Numerical illustrations and discussion 

For illustration of our theoretical results, we have performed numerical so- 
lutions of the Kadanoff-Baym equations including initial correlations. We 
considered the relaxation of electrons in dense hydrogen with self-energies in 
second Born approximation. Starting from an initial nonequilibrium distribu- 
tion, we compared the relaxation for two cases: I) without initial correlations, 
C(t ) = and II), with nonzero initial correlations which were chosen in the 
form of the Debye pair correlation function 

C(q, Pl ,p 2 , t a ) = -YzjSlfWffaXl - f( Pl + q )][l - - q )], (29) 
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Figure 1: Time evolution of kinetic, potential, and total energy for zero (solid lines) and 
nonzero (dashed lines) initial correlations. The initial distribution is an uncorrelated equi- 
librium (Fermi) distribution with T = 10000A' and n = 10 21 cm -3 . 



where / = /(to)- As expected, the presence of initial correlations turns out to 
be important on short times. This becomes particularly clear from analyzing 
the time evolution of the mean potential and kinetic energy, Fig. 1. While 
for the uncorrelated initial state, potential energy starts with zero and builds 
up continuously, the picture changes if there exist initial correlations. With 
the choice of the form (p9[), the correlations are stronger than in equilibrium 
which corresponds to a larger magnitude of potential energy at t = to, which, 
consequently, is reduced in the course of the relaxation. Due to conservation of 
total energy, kinetic energy shows exactly the opposite trend. One clearly sees 
the decay of the initial correlation term, as the curves for the two cases merge 
after times of the order of the correlation time. This confirms that indeed, 
Bogolyubov's weakening principle is reproduced by the presented generalized 
Kadanoff-Baym equations in a dynamic and selfconsistent way. 

Our numerical results illustrate the effect of initial correlations on the 
short-time relaxation behavior for a simple model case. But our theoretical 
approach is completely general and allows for numerical investigations of far 
more complex initial correlations. Besides the fundamental interest in the 
problem of initial correlations in the Kadanoff-Baym equations, our results 
are also of practical importance. With the possibility to start quantum kinetic 
calculations from a general initial state, the scope of nonequilibrium processes 
in many-body systems which are accessible for numerical investigation is es- 
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sentially extended. Although the determination of C(t ) can be complicated 
by itself, our approach allows to separate this problem (the "generation" of 
the correlated state) from the relaxation dynamics. 
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